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00 ■ Abstract 

Binary symmetry constraints of the A/'-wave interaction equations in 1 + 1 and 2 + 1 
dimensions are proposed to reduce the 7V--wave interaction equations into finite-dimensional 
LiouviUe integrable systems. A ne-w involutive and functionaUy independent system of poly- 
' nomial functions is generated from an arbitrary order square matrix Lax operator and used 

to sho-w the LiouviUe integrabihty of the constrained flo-ws of the J\f-wave interaction equa- 
tions. The constraints on the potentials resulting from the symmetry constraints give rise 
to involutive solutions to the TV- -wave interaction equations, and thus the integrabihty by 
quadratures are sho-wn for the 7V--wave interaction equations by the constrained flo-ws. 

' Running title: Symmetry Constraints of Af-wave Equations 

O 

o . 

^ ■ 1 Introduction 
O 

It is a usual practice to utilize the idea of linearization in analyzing nonlinear differential or 
^ ' differential-difference equations (see for example f^, 0] ) • The method of inverse scattering trans- 

_ ^ ■ form is an important application of such an idea to the theory of soliton equations ^] , -which 

^ ' has been recognized as one of the most significant contributions in the field of applied mathemat- 

■ ics in the second half of the last century. The general formulation of Lax pairs is a spectacular 

tool of realization of inverse scattering transform by -which one can break a nonlinear prob- 
lem into a couple of linear problems and then handle the resulting linear problems to solve the 
nonlinear problem. 

Recently in the past decade, an unusual -way of using the nonlinearization technique arose 



in the theory of soliton equations [^-[10|. Although using the idea of nonlinearization is not 
normally considered to be a good direction in studying nonlinear equations, one gradually re- 
alizes that the nonlinearization technique provides a po-werful approach for analyzing soliton 
equations, especially for sho-wing the integrabihty by quadratures for soliton equations. The ma- 
nipulation of nonlinearization not only leads to finite-dimensional LiouviUe integrable systems 



-[15|, but also decomposes infinite-dimensional soliton equations, in -whatever dimensions, 
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into finite-dimensional Liouville integrable systems [16|-|18]- Moreover, it narrows the gap be- 



tween infinite-dimensional soliton equations and finite-dimensional Liouville integrable systems 



1 11, 16, and paves a method of separation of variables for soliton equations |1£, 2C], which 



can also be used to analyze the resulting finite-dimensional integrable systems [^]-[23]. Math- 
ematically speaking, much excitement in the study of nonlinearization comes from a kind of 



specific symmetry constraints [^-|27], engendered from the variational derivative of the spec- 



tral parameter \2t, 27 1. It is due to symmetry constraints that the nonlinearization technique is 



so powerful in showing the integrability by quadratures for soliton equations 29 1. The study 
of symmetry constraints itself is an important part of the kernel of the mathematical theory of 
nonlinearization, which is also a common conceptional umbrella under which one can manipulate 
both mono-nonlinearization Q and binary nonlinearization |26]. 



However, all examples of application of the nonlinearization technique, discussed so far, are 
related to lower-order matrix (here, and in what follows, a matrix is assumed to be square) 
spectral problems of soliton equations, most of which are only concerned with second-order 
traceless matrix spectral problems. On the one hand, there appears much difficulty in handling 
the Liouville integrability of the so-called constrained fiows generated from spectral prob- 
lems, in the case of the third-order and fourth-order matrix spectral problems |2^, ^ It 
is a challenging task to extend the theory of nonlinearization to the case of higher-order ma- 
trix spectral problems. On the other hand, one also notices that mono-nonlinearization can not 
be carried out in the cases of odd-order matrix spectral problems and even-order, including 
the simplest second-order, non-traceless matrix spectral problems. Even for even-order traceless 
matrix spectral problems, it is not clear how to determine pairs of canonical variables to obtain 
Hamiltonian structures of the constrained fiows while doing mono-nonlinerization. Therefore, 
one has to take into account adjoint spectral problems and manipulate binary nonlinearization 
for the case of general matrix spectral problems. In the theory of binary nonlinearization 
there exists a natural way for determining symplectic structures to exhibit Hamiltonian forms 
of the constrained flows. 

In this paper, we would like to establish a concrete example to apply the nonlinearization 
technique to the case of higher-order matrix spectral problems, by manipulating binary nonlin- 
earization for arbitrary-order matrix spectral problems associated with the AA-wave interaction 
equations in both 1-1-1 and 2-1-1 dimensions. The resulting theory will show a direct way for 
generating sufficiently many integrals of motion, and more importantly for proving the func- 
tional independence of the required integrals of motion, for the Liouville integrability of the 
constrained flows resulting from higher-order matrix spectral problems. 



Let us recall some basic notation on binary nonlinearization (see, for example, [33| for a 
detailed description). Let us assume that we have a matrix spectral problem 

= U(p = U{u,X)(j), U = {Uij)rxr, (p = {(pl, - ■ ■ ,fprV (1-1) 

with a spectral parameter A and a potential u = {ui, • • • ,Uq)'^ . Suppose that the compatability 
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conditions 

Ut^ - y^"^^ + ^^""^l = 0, m > 0, 
of the spectral problem and the associated spectral problems 

ct>t^=V^'^^ = V^'^Hu,u,,.---\)^, yM = m>0, (1.2) 

determine an isospectral (At,„ = 0) soliton hierarchy 

5H 

ut^ = X„i{u) = JGm = J^^, m>0, (1.3) 

ou 

where J is a Hamiltonian operator and Hm are Hamiltonian functionals. Obviously, the com- 
patability conditions of the adjoint spectral problem 

V', = -C/^(n,A)V', ij = ii;i,---,Af, (1-4) 

and the adjoint associated spectral problems 

V't™ = -y('")^A = -V^^^^iu, u,,---- A)V (1.5) 

still give rise to the same hierarchy ut^ = X^iu) defined by (pT^). It has been pointed out 
S\ 

|26, O] that J — is a common symmetry of all equations in the hierarchy (1.3). Introducing N 
— ' ' ou 

distinct eigenvalues Xi, X2, • " i ^N, we have 

# = U{u,Xs)<p'^'\ 4'^ = -U''{n,Xs)4'\ l<s<N, (1.6) 

and 

4:l = V^^Hu,u,,---;\s)^^'\ V'2 = -y("^)^(n,n,,---;A,)V'(^), 1 < s < iV, (1.7) 

where we set the corresponding eigenfunctions and adjoint eigenfunctions as (f)^^^ and 1 < 
s < A^. It is assumed that the conserved covariant Gmo does not depend on any derivative of u 
with respect to x, and thus the so-called general binary Bargmann symmetry constraint reads 
as 

Xmo =J2e,^s/-^, i.e., JG^o = jf2f^s¥'^^^^^^^cP^'\ (1.8) 

s=l s=l 

where /x<j, ^ < s < N, are arbitrary nonzero constants, and Es, 1 < s < N, are normalized 
constants. The right-hand side of the symmetry constraint (^]^) is a linear combination of 
symmetries 

ou ou 

Such symmetries are not Lie point, contact or Lie-Backlund symmetries, since i;^^*'^ and 4'^^ can 
not be expressed in terms of x, u and derivatives of u with respect to x to some finite order. 
Suppose that ( |L^ ) has an inverse function 

u = u = 2(</.(^) , . • • , </)(^) ; r--,^''^), (1.9) 
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Replacing u with u in the system (L6) or the system ( |l.7D , we obtain the so-called spatial 
constrained flow: 

# = C/(5,A,),/.W, Vi') = A,)^(^), 1 < ^ < iV, (1.10) 

or the so-called temporal constrained flows: 

^2=yM(s,S,,...;A,)c/>(^), ^t^ = -y(-)^(5,n,,...;A,)V'(^), 1 < . < iV. (1.11) 

The main problem of nonlinearization is to show that the spatial constrained flow ( |1.10| ) and 
the temporal constrained flows ( [1.11 ) under the control of ( |1.10| ) are Liouville integrable. Then 



if (j)^^'^ and 1 < s < A^, solve two constrained flows (|1.10| ) and ( |1.11| ) simultaneously, 

u = u will give rise to a solution to the mth soliton equation ut^ = Xmiu)- It also follows 
that the soliton equation ut^ = X^iu) is decomposed into two flnite-dimensional Liouville 
integrable systems, and u = u presents a Backlund transformation between inflnite-dimensional 
soliton equations and finite-dimensional Liouville integrable systems. More generally, if a soliton 
equation is associated with a set of spectral problems 

0,, = [/«(n,A)</>, l<i<p, 

then it will be decomposed into p + 1 flnite-dimensional Liouville integrable systems. The above 



whole process is called binary nonlinearization |16, 33 1. 



This paper is structured as follows. In Section 2, we will present binary symmetry constraints 
of the A/'-wave interaction equations in 1 -|- 1 dimensions, and show Hamiltonian structures and 
Lax presentations of the corresponding constrained flows. In Section 3, we consider the 2-1-1 
dimensional case. We will similarly construct binary symmetry constraints of the A^-wave interac- 
tion equations in 2-1-1 dimensions, and discuss some properties of the corresponding constrained 
flows. In Section 4, we go on to propose an involutive system of functionally independent polyno- 
mial functions, generated from an arbitrary-order matrix Lax operator, along with an alternative 
involutive and functionally independent system. An r-matrix formulation will be established for 
the Lax operator, and used to show the involutivity of the obtained system of polynomial func- 
tions, together with Newton's identities on elementary symmetric polynomials. A detailed proof 
will also be made for the functional independence of the system of polynomial functions by using 
the determinant property of the tensor product of matrices. In Section 5, two applications of 
the involutive system engendered in Section 4 will be given, which verify that all constrained 
flows associated with the A/'-wave interaction equations in both 1-1-1 and 2-1-1 dimensions are 
Liouville integrable. Moreover, a kind of involutive solutions of the AA-wave interaction equations 
in two cases will be depicted. Theses also show the integrability by quadratures for the A^-wave 
interaction equations. Finally in Section 6, some concluding remarks will be given, together with 
conclusions. 



2 Binary symmetry constraints in 1 + 1 dimensions 
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2.1 nxn AKNS hierarchy and 1 + 1 dimensional A/^-wave interaction equations 

Let n be an arbitrary natural number strictly greater than two. We begin with the nxn matrix 
AKNS spectral problem ||] 

0,. = [/(/> = U{u, X)^, U{u, X) = XUo + = (01, • • • , ^nf, (2.1) 

with a spectral parameter A and 

[/o = diag(ai, • • • ,a„), C/i(n) = (njj)„xn, (2.2) 

where Oj, 1 < i < n, are distinct constants, and uu = 0, 1 < i < n. The standard AKNS 
spectral problem, i.e., the spectral problem (p.l|) with n = 2, has been analyzed in |35|, but it 



can not generate any A/'-wave interaction equations and thus it is not discussed here. In order 
to express related soliton equations in a compact form, we write down the potential u as 

J -u = p([/), i.e., n = (ii2i,ni2,'U.i3,'"3i,''i23,'"32)^, when n = 3, 

I (2.3) 
,ui2,ui^,u:ii,ui4,U4i,U2s-,u^2i ' ' ' i'Un,n—ii'iJ-n~i,n) i when n > 4, 

in which we arrange the exponents Uij in a specific way, first from smaller to larger of the integers 
k = i+ j and then symmetrically for each set < i < A; — 1}. 

Let us now consider the construction of the 1 + 1 dimensional AA-wave interaction equations 
and its whole isospectral hierarchy associated with the spectral problem (|2.1| ). We first solve the 
stationary zero-curvature equation for W: 

- [U, W]=0, W= {W^j)nxn, (2.4) 

which is equivalent to 

n 

ll^ij> + Uij{Wii - Wjj) + ^ {ukjWik - UikWkj) - X{ai - aj)Wij = 0, i ^ j, 



k=i 



(2.5) 



Wii^x = "^{uikWki - UkiWik) 



fc=i 

where 1 < i,j < n. We look for a formal solution of the form 

.(0 

and thus (|2.5|) becomes the following recursion relation 



W = Y, WlX-', Wl = {Wif )ny<n, (2.6) 



n 

k = l 

n 

k = l 



(2.7) 



where 1 < i,j < n and / > 0. In particular, from the above recursion relation, we have that 



W^f =Pi = const., = 0, 1 < i / i < n, (2i 



and 



We require that 



W^' = 0, Wll' = ' \_ ' uij, 1 < i / j < n. (2.C 



w'>f\u=o = ^, l<i,j<n, />!. (2.10) 



This condition (^lo| ) means to identify all constants of integration to be zero while using (p. 
to determine VF, and thus all Wi, I > 1, will be uniquely determined. For example, we can 
obtain from ( |2.7D under ( |2.10| ) that 



Wy = —Uij^x H > , ( '-)UikUkj, I J <n 



[Oi - ajY Oi - Oj ^ Ok - Oi Ok - aj 

(2.11) 

Ti^(2) sr^ Pk - Pi 1 / • / 

fc — 1 



It is easy to see that the recursion relation (2/7) can lead to 



k = l 

n 

+ [u^J^-WjWf^ - {uik + Uijd-^ujk)wj!^^] = (a, - a,)W^.^^\ i + j, 



1 



(2.12) 

where 1 < i,i < n, ^ > 1, and d^^ is the inverse operator of 5 = This can be written as the 



Lenard form 

MG/_i = JGu I > 1, (2.13) 

where Gi = p{Wi+i) is generated from M^^+i in the same way as that for u, and J is a constant 
operator 



J = diag(^(ai - 02)0-0, (ai - 03)0-0, (02 - a3)o-oj , when n = 3, 

J = diag(' (ai - 02)0-0, (ai - 03)0-0, (ai - 04)0-0, (02 - 03)0-0, • ■ ■ , («n-i - Qn)o-oJ , (2.14) 

n(n-l)/2 

when n > 4, 
with 0-0 being given by 

-(:,;) 
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For example, when n > 4, we have 

= {WSMI wS,wflwS,wfl W^l W!il . . . , W^i'li.n)^, I > 1, (2.15) 

the first of which reads as 

//3i-/32 /3i-/32 /?i-/33 /?i-/33 /3i-/34 /3i - /?4 

(jo = ^^21, ■U12, ^il3, ■"41, WM, 

Vai — 02 ai — 02 «! — Q!3 «! — as ai — a\ — 

-Un,n-1, Un-l,n] ■ (^-loj 



Q^n— 1 Ol-n OL-n—l Ot% 

The operators J and M are skew-symmetric and can be shown to be a Hamiltonian pair 

We proceed to introduce the associated spectral problems with the spectral problem ( ^ ) 

= T/M</,, T/M = T/M(^x,A) = (A'"Ty)+, m > 1, (2.17) 

where the symbol + stands for the choice of the part of non-negative powers of A. Note that we 
have 

Wi^ = [Uo,Wi+i] + [U^,Wi], />0, 

and we can compute that 

m 
1=0 

m m 

1=0 1=0 

m—1 m 

= ^[C/o,VFm]A™-' + ^[C/i,VF^]A™-', 

1=0 1=0 

where we have used [C/q, Wq] = 0- Therefore, under the isospectral conditions 

At™. = 0, m > 1, (2.18) 
the compatibility conditions of the spectral problem ( |2.1[) and the associated spectral problems 



( 2.17 ), i.e., the zero-curvature equations 
equivalently lead to 

Uu^ = Wrax - Vx.Wm\ = Wo, VF^+i], m > 1. 

This gives rise to the so-called n x n AKNS soliton hierarchy 

ut^ = Xm ■■= JGn,, m>l, (2.19) 
where J and Gm = piWm+i) are determined by ( 2.14| ) and ( |2.13| ). 
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Applying the trace identity 

6u 

where 7 is a constant to be determined, we can obtain 



Sua 



1 



(2.20) 



in which 1 < i ^ j < n and 7 is determined to be zero. In this computation, we need to note 
that 

tr(H^— ) = tiiWUo) = Y^iaiWS + a2wS + ■■■ + a„W«)A"^ 



l>0 



and 



dU 



tr(Vr^) = tv{WEij) = Wji = ^ W^lh~\ 1 < i / J < n. 



l>0 



where Eij is an n x n matrix whose entry is one but other entries are all zero. Therefore, 
the isospectral hierarchy ( 2.19 ) has a bi-Hamiltonian formulation 



Xr, 



J- 



M— — , m > 1. 



(2.21) 



6u 6u 

The first nonlinear system in the hierarchy ( ^.19 ) is the 1 + 1 dimensional A/"- wave interaction 



equations [39| 



Pi - Pj , sr^ , Pi - Pk Pk - Pj 1 / ■ / • / 
-Uij^x + > , ( -jUikUkj, I J <n. 



This system is actually equivalent to the following equation in the matrix form 

Uit, =Wi,-[Ui,Wi], 

which can be rewritten as 

Pt,=Q.-[P,Q], [Uo,Q] = [Wo,P], 



(2.22) 



(2.23) 



(2.24) 



where P and Q are assumed to be two off-diagonal potential matrices. Based on ( 2.231) , a vector 
field p{SP) is a symmetry of ( |2.22| ) if the matrix 6P satisfies the linearized system of ( ^.22]) : 



{6P)t, =i6Q),-[Ui,dQ]-[6P,Wi 



with SQ being determined by 



[Uo,5Q] = [Wo,SP]. 



(2.25) 



(2.26) 



The AA-wave interaction equations ( |2.22D contains a couple of physically important nonlinear 
models as special reductions |^^, for example, three- wave interaction equations arising in fluid 
dynamics and plasma physics |41, 42, with U being chosen to be an anti-Hermitian matrix. 
Its Darboux transformation has been established in [Q], which allows one to construct soliton 
solutions in a purely algebraic way. The Darboux transformation has also been analyzed for the 
AA-wave interaction equations with additional linear terms [^]. 
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2.2 Binary symmetry constraints in 1 + 1 dimensional case 



We would like to present binary symmetry constraints of the 1+1 dimensional A/"- wave interaction 
equations (|2.22| ). To this end, we need to introduce the adjoint spectral problem of (^H]): 



(2.27) 



and the adjoint associated spectral problem of (|2.17| ): 

^pt^ = -V'^'^^''{u,X)^P, (2.28) 

where U and are given as in ( p.l| ) and (2.17), respectively. The compatability condition of 

(|2]2^) and (|]2|) still gives rise to ut^ = defined by (^). 

The variational derivative of the spectral parameter A with respect to the potential u can 
be calculated by (see ||2^, 28 1, or [16 1 for a detailed deduction) 



5\ 
6u 



du 



I.e., 



5\ 
Sua 



E 



-1, 



iijjj, 1 < i / j < n, 



(2.29) 



where E is the normalized constant: 



ip -Trr(pdx. 
3 OA 

A direct calculation can show that the variational derivative satisfies the following equation 



M 



5\ 
6u 



XJ 



6X 

Su 



(2.30) 



Since A does not vary with respect to time, we have a specific common symmetry of the 



hierarchy ( 2.19| ). To carry out binary nonlinearization, we take a Lie point symmetry of the 
AA-wave interaction equations ( |2.22| ), 



Yo:= pi[T,Ui]), r = diag(7i,. 



.7nJ 



(2.31) 



where 71,72, ■'■ )7n are arbitrary distinct constants (Xq = JGq is an example with T = Wq). It 
can be easily checked that 

{6P,6Q) = {[r,Ui],[r,Wi]) 

satisfies ( |2.25| ), and thus Yq is a, symmetry of ( p.22D . Then, make the following binary Bargmann 
symmetry constraint 



Yo = fiEJ-- = ^JV;^ 

ou ou 



dU 



(2.32) 



where fj, is an arbitrary nonzero constant, J is defined by ( 2.14| ), and (f) and ip are the eigen- 
function and adjoint eigenfunction of (2.1) and (2.27), respectively. Upon introducing N distinct 
eigenvalues Ai, A2, • • • , Aat, we obtain a general binary symmetry constraint 



(2.33) 
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where fis, ^ ^ s < N , are N nonzero constants, and (f>^^^ and 1 < s < A^, are eigenfunctions 
and adjoint eigenfunctions defined by 



and 



= y(i)(n,A,)</)W, 4f = -y(i)^(n,A,)V^(^), 1 < s < iV. 
Let us rewrite the left-hand side of ( 2.33| ) as the matrix form 

N 

6P = p-\Zo) = [Uo,J2f^s^^'^i'^'^''h 



(2.34) 
(2.35) 



(2.36) 



s=l 



which allow us to prove, by a direct computation as in |4(:] but more conveniently, that the 
vector field Zq = p{5P) is really a symmetry of the A^-wave interaction equations ( p.22 ). Now 
the symmetry problem is equivalent to showing that 



N 



N 



(2.37) 



satisfies the linearized system (2.25), when (f)^^^ and V'''*'') ^ ^ s < N , satisfy (2.34) and ( |2.35| ). 
A detailed proof will be given in Appendix A. 

Therefore, we have the following binary symmetry constraint 

Yo = jf;^.V'(^)^^^^^0(^\ i.e., [r,C/i] = [Uo,J2f^s^^'^4'n (2.38) 



s=Q 



s=l 



When N and ps vary, ( |2.3^ ) provides us with a set of binary symmetry constraints of the A/"- wave 
interaction equations (|2.22|). Let us assume that 



(2.39) 



= {cpis, (t>2s, (t>nsf , tp^"^ = {tpls,tp2s, • • • , tpns)'^ 



in order to get an explicit expression for u from the symmetry constraint ( 2.38| ), and introduce 
two diagonal matrices 



A = diag(Ai, ■ . . , \n), B = diag(//i, . • ■ ,iin) 



(2.40) 



which will be used throughout our discussion. Solving the Bargmann symmetry constraint ( |2.38| ) 
for n, we obtain 



Q,- — 



li - Ij 



1 j <n, 



where B is given by ( 2.40| ), and and are defined by 

$i = ((/)ii,(/)i2, • • • ,(/'j7v)^, ^'j = (■i/'ii,'0j2, • • • jV'jiv)^, l<i<n, 
and (., •) denotes the standard inner-product of the Euclidean space R^. 



(2.41) 



(2.42) 
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Note that the compatabihty condition of ( 2.34 ) and ( 2.35 ) is still nothing but the 1 + 1 
dimensional A^-wave interaction equations ( ^.22 ). Now using ( 2.41 ), we nonlinearize the spatial 
part ( |2.34| ) and the temporal part ( 2.35| ) of spectral problems and adjoint spectral problems of 
the AA-wave interaction equations ( 2.22| ). Namely we replace Uij with Uij in N replicas of the 



spectral problems and adjoint spectral problems ( 2.34 ) and N replicas of the associated spectral 
problems and adjoint associated spectral problems ( |2.35| ), and then obtain two constrained flows 



for the AA-wave interaction equations (2.22): 

# = Uiu,Xs)^^'\ = -U^iu,Xs)4'\ l<s<N, 



and 



where u = p{{uij)nxn) is defined like u. For example, when n > 4, we have 



(2.43) 



(2.44) 



(2.45) 



In order to analyze the Liouville integrability of the above two constrained flows, let us first 
introduce a symplectic structure 



n N 



(2.46) 



i=l i=l s=l 

over M^"^, and then the corresponding Poisson bracket 



i=l 



n N 



dipis d4>is dcpis dipi. 



(2.47) 



i=l s=l 

where the vector field Idf is defined by 

u\X,Idf ) = df{X), X G r(M2«^). 

A Hamiltonian system with a Hamiltonian H defined over the symplectic manifold (R^"'^,a;^ 
is given by 



= {^^,H} = -B-'^, = {^i,H} = B-' — , l<i<n, 



(2.48) 



where t is assumed to be the evolution variable. Second, we need a matrix Lax operator 

L«(A) = Ci + Di(A), (2.49) 

with Ci and Di{\) being defined by 



Ci = r = diag(7i,---,7n), I)i(A) = (Z);;^(A))„xn, I)g^(A) = ^ 



N 



A- A 



-(k^si^js, (2.50) 
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where 1 < i,j < n. Note that upon taking binary nonlinearization, we obtain 

a,- — a,- 



U{u,X) = XUo + Ui{u) = XUo + (uij), Uij 
V^^\u, A) = A^o + Wi{u) = XWo + (vij) 
where 1 < i,j < n. 



li - Ij 

._ (3i - I3j ^ _ (3i- 



(2.51) 



[^^,B^J), (2.52) 



Theorem 2.1 Under the symplectic structure (2.4('), the spatial constrained flow (2.4j) and 
the temporal constrained flow ( 2.44 ) /^'^ 1 + 1 dimensional Af-wave interaction equations 
( 2.22J are Hamiltonian systems with the evolution variables x and ti, and the Hamiltonians 



k=l 



l<k<l<n 



k=l 



l<k<l<n 



ak - ai 
Ik - li 

Pk-Pi 
Ik -li 



{^k,B^i){<^uB^k). 



(2.53) 
(2.54) 



respectively, where A and B are defined by ( ^.4C\) , and and ^i, 1 < i < n, are defined by 
( 2.4^) - Moreover, they possess necessary Lax representations, i.e., we have 



(lW(A)), = [t/(5,A),L«(A)], (L«(A))t, = [y«(5, A), L«(A)], 



(2.55) 



where L'^^^X), U, and V'^^\X) are given by ^4^, \2.5(\ ), and ( \2.5^ ), if and ^2.4A ) 

hold, respectively. 



Proof: A direct calculation can show the Hamiltonian structures of the spatial constrained 
flow ( 2.43| ) and the temporal constrained flow (2.44) with Hf and defined by ( p. 53 ) and 



( 2.54 ). Let us then check the Lax representations. By using ( 2.43 ), we can compute that 

N 



s=l 



N 

E {U{u, A.)</.W V^W^ - <^W^W^C/(S, A.)) 



s=l 

N 



s=l 



A- A 



TV 



[C/(a,A),L«(A) - Ci] - [C/o,^/^^'/'^'^^^^'^] 

s=l 

N 

[U{u,X),lW(X)] + [Ci,U{u,X)] - [Uo,J2l^scl>^'^^^'^'^] 

s=l 

N 

[[/(S,A),L«(A)] + [Ci,Ui{u)] - [Uo,Y.f,s<f>^'^i^^''^^]. 



s=l 
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This implies that (L(i)(A))^ = [U{u, X), L^^^X)] if and only if 



N 



s=l 

The above equality equivalently requires the constraints on the potentials shown in ( |2.41| ). 
Therefore, the spatial constrained flow ( 2.43| ) has the necessary Lax representation defined as in 



(|2.55| ). The proof of the other necessary Lax representation (^(-'^^(A))^^ = [y*^^''(n, A), ^(-'^^(A)] is 
completely similar, and thus we omit it. The proof is finished. I 



We remark that the Lax representations ( 2.55 ) are not sufficient. Namely, we can not obtain 



the spatial constrained flow ( |2.43| ) or the temporal constrained flow ( 2.44 ) from the correspond- 



ing Lax representation in ( p. 55 ). This can be easily observed by considering a special class of 



solutions of ( |2.55| ). For example, either any vector functions i;^^''^ with ■^^^^ = 0, 1 < s < A^, or 
any vector functions ip^'^^ with (p^^^ = 0, 1 < s < N , will solve ( 2.55 ), but it is easy to see that 



they do not always solve (|2.43D [or (|2^ )] since (f)^^^ and 1 < s < N, have to solve some 
ODEs resulting from ( |2l^ ) [or {^A^)]. 



3 Binary symmetry constraints in 2 + 1 dimensions 
3.1 2 + 1 dimensional AT- wave interaction equations 

Let n be an arbitrary natural number strictly greater than two. Similar to the case of the 1+1 
dimensional A/'-wave interaction equations, let us begin with the Lax system 

Fy = JF, + PF, Ft = KF, + QF, F = (/i, • ■ ■ , (3.1) 

in 2 + 1 dimensions. Here it is assumed that 

J = diag(Ji, • • • , J„), K = diag(Eri, • • ■,Kn), Ji / Jj, Kij^Kj, 1 < i / j < n (3.2) 

are two constant diagonal matrices, and P and Q are two n x n off-diagonal potential matrices 

P = P{x, y, t) = {Pij)nxn, Q = Q{X, y, t) = {qij)nxn- (3.3) 

The compatability condition Fyt = Fty of the Lax system ( |3.lD reads as 

[J, Q] = [K, P], Pt-Qy + [P, Q] + JQ, -KP, = 0, (3.4) 

which is called the 2 + 1 dimensional J\f-wave interaction equations [^ . The equation [J, Q] = 
[K, P] tells us that Q can be represented by P and vice versa, and so practically, we have just 
one of two potential matrices to be solved. The adjoint system of the Lax system ( |3.1[ ) is given 
by 

Gy = JG, - P'^G, Gt = KG, - Q^G, G = (51 , ■ • • , Qnf , (3.5) 
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whose compatability condition Gyt = Gty still gives rise to the 2 + 1 dimensional TV-wave 



interaction equations (|3.4[) . 

We first use a symmetry constraint of the 2 + 1 dimensional A/'-wave interaction equations 
to change the above problem in 2 + 1 dimensions to three problems in 1 + 1 dimensions. 
], we introduce the spectral problems 



As made in [48 



: G, X)(l) = (Ar?g + nf)(P = 

■ ny{p, F, G, A)(^ = {xnl + n\)(t> 
n\Q, F, G, A)(/> = {xnl + - 



A/„ F 
G^ 



XJ + P JF 



G^J 







(3.6) 



XK + Q KF 
G^K 



where In is the nth-order identity matrix and = ■ • • , (pn+i)'^ ■ The new extended po- 
tentials in the above spectral systems consist of not only the original potentials, P and Q, but 
also the solutions of the Lax system and the adjoint Lax system, F and G. The compatability 
conditions cpxy = (t'yx-, <Pxt = <Ptx, and (pyt = (j)yt give rise to the 2 + 1 dimensional TV- wave 
interaction equations (3.4), the original Lax system ( p.lj ) and its adjoint system (3.5), and the 
nonlinear symmetry constraint of (|3.4D: 



= [FG^,J], Q, = [FG' ,K]. 



T 



(3.7) 



It is easy to check that {6P,6Q) = {[FG'^, J], [FG'^ , K]) satisfies the linearized system of the 
2 + 1 dimensional TV- wave interaction equations ( 0.41) : 

[J, 6Q] = [K, 5P], {5P)t - {5Q)y + [5P, Q] + [P, 6Q] + J{6Q), - K{5P), = 0, (3.8) 

when F and G solve the Lax system ( |3.lD and the adjoint Lax system ( |3.5| ), respectively. There- 
fore, (|3.7D is really a symmetry constraint of the 2 + 1 dimensional TV- wave interaction equations 
(|3.4| ), since both sides of ( |3.7D are symmetries of (p.4|). Now we see that the original problem in 
2 + 1 dimensions is transformed into three problems in 1 + 1 dimensions. The spectral problems 
are our starting point to make a link of the 2 + 1 dimensional TV- wave interaction equations 
to finite-dimensional integrable systems. 



3.2 Binary symmetry constraints in 2 + 1 dimensional case 

Let us start from the spectral problems in (^]^), which are similar to those for the 1 + 1 dimen- 
sional TV- wave interaction equations ( p. 22]) . The main difference is that the coefficient matrix of 
A in the x-part of the spectral problems (|3.6| ) is 

f)g = diag(l,---,l,0), (3.9) 
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whose diagonal entries are not distinct. However, the y-part of the spectral problems ( |3.£| ) has 
the same property as the spectral problem (|2.l| ) in 1 + 1 dimensions. Therefore, we use the y-part 
of the spectral problems (l3.6|) to compute the variational derivatives of A: 



6pij dpij 6qij dqij Ki - Kj 

— = E-'ij^--^ = E-'j^cPn+iA, J- = E-^^^^^ = E-'Jict)i^n+i, l<i<n, 

where E is the normalized constant, and tp = (V'l, • • • , V'nj is an adjoint eigenfunction of 

the adjoint spectral problems 

V-. = G, A))^V = -im^of + i^ff)i^ = -[fT ) V', 

i;y = -{nviP, F, G, A))V^ = -iXinl) + inyf)i; = - ( ^-^Jrf "^(f ) V', (3.10) 

= -{n\Q, F, G, A))^V = -imif + {^{f)^ = - ( ) 

These variational derivatives of A give us a conserved covariant and also a clue to compute a 
required symmetry, expressed in terms of eigenfunctions and adjoint eigenfunctions. 

As in the 1 + 1 dimensional case, upon introducing N distinct eigenvalues Ai, A2, • • • , Atv, 
we have 

4^) = r!"(n,A,)(/.(^), </)W = ny{u,Xs)4>^'\ 4>f = n\u,Xs)(t>^'\ i<s<n, (3.11) 

and 

= -(i^^-f (n,A,)v^w, = -{nyf{u,Xs)^^'\ 4'^ = -{nY{u,Xs)4'\ i<s<n, 

(3.12) 

where i;^^'^^ and n + 1 dimensional vector functions: 

= ((/.I,, • • • , (^ns, (^n+l,sf, = ii^ls, • • • , Vns, V'n+l,s)^, 1 < S < N. (3.13) 

To carry out binary nonlinearization, we need to construct two special symmetries, the one of 
which is a Lie point symmetry, and the other of which is not a Lie point, contact or Lie Backlund 
symmetry, but generated from ( p. 11 ) and (3.12). Let us choose a set of n + 1 arbitrary distinct 
constants 6i, - • • ,6n, Sn+i, and set 

A = diag(5i,---,5„). (3.14) 
Similar to the 1 + 1 dimensional case, it can be directly shown that 

i6P,5Q,6F,dG) = i[A,P],[A,QiAF - 5n+iF,AG - 5n+iG) (3.15) 

and 

f Spij = (J, - Jj){<l>i,B^j), 6qij = {Ki - Kj){^u B^,), 1 < i / i < n, 

< (3.16) 

I 6f^ = {<l>i,B^ri+l), Sgi = {<^n+l, B^i), 1 < 1 < U, 
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are two symmetries of the equations (|3.4|) , (13. if ) and (|3.5D . That is to say that they satisfy the 
hnearized system of the equations (0) , (|3.l|) and (|3.5|) : the first subsystem (|3.8D and the second 
subsystem 



i6F)y = J{6F), + i5P)F + P5F, {5F)t = K{6F), + {5Q)F + Q6F, 
{6G)y = J{6G)^ - {6PfG - P'^5G, {6G)t = K{5G)^ - {dQfG - Q^SG, 

for all solutions {P,Q,F,G) of Q, and Q. Here we remind that 



(3.17) 



B = diag(/ii, • • ■ ,fiN)'^ 



is defined by ( p.40 ), (•, •) denotes the standard inner product of ffi'^, and $j and are similarly 
defined as 

= {(j)ii,(pi2,--- ,4>iN)'^, "^i = {ipii,'>Pi2, - ■ ■ , l<i<n + l. (3.18) 

Now a binary Bargmann symmetry constraint of (|3.4|), (|3.1|) and (p.5|) can be taken as 



([A,P]),,- = (J,- J,)(cl>„i?^j), ([A,Q]),, = iKi-K,){<l>i,B^,), l<i^j<n, (3.19) 
(AF - 6n+iF)i = (AG - 6n+iG)i = 51-^), 1 < i < n. (3.20) 

This symmetry constraint gives us the following choice for the constraints on the extended 
potentials 

Pij=Pij ■■= ^^^^{^^,B^^>,), qi, = q,, := ~ {'^i, B^j), I < i ^ j < n, (3.21) 

Oi Oj Oi Oj 

h = h-=^-k ($i,5*n+i), 5* = 5* := ($n+i,5**), l<i<n. (3.22) 

Oi — On+l Oi — On+1 

One can express the above symmetry constraint in another way. Actually, it can be proved that 

(<5P,,5g) = ([A,P],[A,Q]), 

and under the constraint (3.22), 



5Pij = {Ji - Jj){^^, Sqij = {K, - Kj){^„ B^j), l<i^j<n, 

are two symmetries of the 2 + 1 dimensional A/"- wave interaction equations (|3.4|). 



Now plug the above expressions for the extended potentials, ( p. 21 ) and ( |3.22| ), into the 



spectral problems ( |3.6| ) and the adjoint spectral problems ( |3.10| ), and then we get the constrained 
flows 

# = n-{F, G, A,)0(^), V'i') = -{n-{F, G, Xs)fi^^'\ (3.23) 
= ny{P, F, G, Xs)4'\ = -{^"{P, F, G, Xs)fi^^'\ (3.24) 
4'^ =n\Q,F,G,Xs)^^'\ # = -(0*(Q,F,G,A,))V), (3.25) 
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where 

P = {Pij)nxn, Q = {qij)nxn, F = {fi, ■ ■ ■ Jnf , G = [qi, ■ ■ ■ ,gnf ■ (3.26) 

All these three constrained flows are systems of ordinary differential equations of (pis and ipis, 
l<i<n + lA<s<N. 



n+l N 



We introduce the symplectic structure 

n+l 

^2 = ^ Bd^i A d^-i = ^ ^ nsdcpis A dtPis (3.27) 

j=l j = l s = l 

over M2(n+i)Af ^pj^g corresponding Poisson bracket and the corresponding Hamiltonian form with 
the Hamiltonian H and the evolution variable t are similarly taken as 

n+l 

{tT fin /IT flfl 

(3.28) 



Similar to Theorem 12. 11, we have 



= m = 5-1^ 1< i < n + 1. 



(3.29) 



Theorem 3.1 Under the symplectic structure (3.21), three constrained flows (3.25), ( 3.2.^ and 
( 3. 23^ ) are Hamiltonian systems with the evolution variables x, y and t, and the Hamiltonians 

n n ^ 
- B^>k) - V ^ ^ {^k, B^>n+l){^n+l,B^>k), (3.30) 



-"2 



l<fe</<n 



{^k,B^i){<^i,B^k) 



Jk 



Sk - 



'n+l 



($fc,5^'„+l)($„+l,S^'fc), 



k=l 

n 

E 

k=l 

Y,I<k{A<^k,B^k)- Yl ^,''~f' {'^k,B^i){'^i,B^k> 



(3.31) 



k=l 

n 

E 

k=l 



l<k<l<n 



Kk 



Sk - 5.. 



'n+l 



(3.32) 



respectively, where A and B are defined by ( 2.4(\ ), and ^i, 1 < i < n + l, are defined by 
( 3.1^ ). Moreover, they possess the necessary Lax representations 

(L(2)(A)), = [17-(i?,G,A),L(2)(A)], (3.33) 
(L(2)(A)), = [17^(P,F,G',A),l(2)(A)], (3.34) 
(L(2)(A)), = [1^*(Q,F,G,A),L(2)(A)], (3.35) 

respectively, where P, Q, F and G are given by ( ^.26i) , (^.21\ ) and l\3.2^) , and L^'^\\) is defined 
by 

{ L(2)(A) = C2 + D2iX), C2 = diag(A, 5n+i) = diag(5i, ■ ■ ■ , 5„, 5„+i), 

^ (3 36) 

D2 = (i^Jf )„+!,„+!, = Y YT^Mos, 1 < i, j < n + 1. 
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Proof: It can be verified by a direct calculation that all three constrained flows ( |3.23| ), 
( 3.24 ) and ( 3.25 ) have the Hamiltonian structures under the symplectic structure ( 3.27| ) with 
the Hamiltonian functions i/f , -f^l ^2 shown in ( ^.30D , ( 3.31D and ( 3.32| ). Let us now check 
three Lax representations ( |3.33 ), ( 3.34 ) and ( |3.35D . Since the proofs are similar for all three 
cases, we just show the second case, i.e., the Lax representation of the constrained flow ( |3.24| ). 
By using (|3.24|) , we can compute that 



N 



A- A 



N 

Y. 

s=l 
N 

E 

s=l 
N 

E 

s=l 



A -A, 

X-Xs 
A -A. 



ny{p, F, G, A,)0Wv^'^^ - (i)'''^'-'^^ny{p, f, g, a,) 



{[ny{p, F, G, A), <a(^)V'«^] - [ny{p, f, g, a) - ny{p, f, g, a,), <a(^)V'(^)^]) 



N 



= [ny{p,F,G,x),L(^\x) - C2] - [i7^,^M.</.(^)^w^] 

s=l 

N 

= [ny{p,F,G,x),L^'\x)] - [nl{p,F,G),C2] - [i^^,^]/..^^^^^^^^]. 

s=l 

Therefore, it follows that {L(^\X))y = [^y {P , F , G , X) , L^-^^X)] if and only if 

TV 

[C2,ny{p,F,G)] = [ny,Y,t^s<P^'^^^''^^]- 



.8 = 1 



This equality equivalently requires the nonlinear constraints on the potentials defined by ( |3.21| ) 
and ( 3.22D . Therefore, the constrained flow ( p.24 ) has the necessary Lax representation shown 
in ( |3.34| ). The proof is finished. I 

We also remark that the Lax representations (|3.33| ), ( 3.34 ) and ( 3.35 ) are not sufficient 
to generate the corresponding constrained flows defined by ( ^.23|) , (^.24 ) and ( 3.25]) , since the 
Gateaux derivative operators of the Lax operators 17^ and given in (|333|) , (1^) and 



j.25| ) are not injective. However, it will be shown that they are good enough in generating 
integrals of motion of the constrained flows. 



4 An involutive and functionally independent system of poly- 
nomial functions 



Let m be an arbitrary natural number. We start from an m-order matrix Lax operator 

L(A) = L(A; ci,---,Cm)=G + D{X), (4.1) 
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with C and D{X) being defined by 

N 



C = diag(ci, • • • ,Cm), D{\) = {Dij{\))mxm, Ai(A) = V'- — ^—(pisipjs, l<i,j<m. (4.2) 

A — As 

s=l 

Here Cj, A^, and Hs are arbitrary constants satisfying 

N 

n^s/O, Xi^Xj, l<iy^j<N, (4.3) 

s=l 

and (/>j<j and il^js are pairs of canonical variables of the symplectic manifold (M^'"^,u;^) with the 
symplectic structure 

m N 

.2 



AdV'is- (4.4) 



i=l s=l 



The corresponding Poisson bracket reads as 

771 N 

4.1 r- matrix formulation 

As usual, two special matrices defined by the tensor product of matrices are chosen as 

Li(A) = L(A) (g) Im, L2{n) = /m «) L{fi), (4.6) 

where is the mth-order identity matrix, and 

{A (g) B)ij^ki = aikbji ifA = (aij) and B = {hj). (4.7) 

We want to find an x matrix r = r(A, /u) so that we have a r-matrix formulation |5C, ^] 

{L(A) ® L(^)} = [r(A,^),Li(A) +L2(m)], (4.8) 

with the Poisson bracket {L{X) f" L{p)} being defined by 

({L(A) ? L(^)}),j,H = {L,k{X),Lji{fi)} = uj\ldLji{fi)JdL,k{X)), 1 < i,j,k,l < m, (4.9) 

where L = {Lij)mxm is assumed. Let us first compute {Ljj(A), Lh(^)}- When i ^ I and j ^ k, 
it is easy to obtain {Lij{X),Lj.i{ii)} = 0. When i ^ I and j = k, we have 



N 



{U,{X),L^i{f,)} = Y, 



^ A — As /i — As 

^ 1 1 

E ^(lA ^H^si^ls = ALaW - LM). 

;U — AA — As fJ- — As fj, — A 
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Similarly, when i = I and j ^ k, we have 

N 



and when i = I and j = k, we have 



N 



{L,, (A), L,,(m)} = E ^ _ ^ _ ^ - A - A, M - A. 



AT 

E 



fi — X 



[{LuiX) - Luiii)) - (L,,(A) - L,-,(m))]. 



Therefore, we obtain 



r 0, 



{Ljj(A),Lfe/(/x)} = < 



when z 7^ Z, j ^ k 

when z = Z, j ^ k 

j;^{Lii{X) - Liiifi)), when z 7^ Z, j = A; 

. 7irA[(^ii(A) - LM) - (LjjiX) - Ljjin))], when i = l, j = k 
In view of this property, we claim that 



(4.10) 



m 



(4.11) 



p,q=l 



where is an m x m matrix with the (p, q) entry being one but the others, zero. Let us second 
compute that 

{[-^V,Li{X) + L2{ix)])ij,ki 

jJi — A 

{P,Li(A)] + [P,L2(m)]W/ 



/X — A 



— E ([E^pg, L(A)] ® Eqp + Eqp [Epq, L{ll)\)ij^kl 



^ — X 



p,q=l 

m 



XI ([-^P<?'-^(A)]ifc(-E'gp)j7 + {Eqp)ik[Epq, L{fi)]jl 



IJ, — X 



H — X 



p,q=l 



{[Eij,L{X)]ik + [Eki,L{n)]ji, 



where we have used {A (g) B){A' (g) S') = (AA') (g) (SB'). Further noting that 

gth 

••• -Lip ■■■ 



[Epq, L] — EpqL — LEpq — pth. 



Lql ' ' ' -^qq Lpp ■ ■ ■ Lq^ 



-L. 



mp 
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we have 



([ rV,Li{X)+L2{fi)])ij,kl 

/X — A 

0, when i ^ I, j k 

-^{Ljk{X) - Ljkifi)), when i = I, j ^ k 

■jj^{-Lii{X) + Lii(fi)), when i ^ I, j = k 

. jrxii^ni^) - Lii{X)) + iUiifi) - Ljjifi))], when i = l, j = k. 
Now ( [4.101 ) and (4.12) shed right on the following theorem. 



(4.12) 



Theorem 4.1 If L(X) = L{X;ci, - ■ ■ ,Cm) is defined by (4-1) and (4-^), then the r- matrix for- 
mulation 



{L{X)^ L{^i)} = [r{X,n),L{X)0lm + Im®Li^i)], r = Eij ® Eji (4.13) 

fJ, — A '-^ 

holds for arbitrary constants ci , C2 , • • • , Cm, • 



It follows from ( 4.13 ) that 

{L\X) ? L\ii)} = [r^'''(A,M),^i(A) +L2(/u)], k,l>l, 
where r'^''(A,^) is given by |52] 



k I 



r''\X,^^) =^^L'-\X)L^i\f,)riX,f^)L\-\X)Ll\^,). 
i=i j=i 

Since for A = {aij)mxm and B = {bij)mxm we have 

m m 

tr {A^B]= ^ {A ? B}ij,i, = bj,] = {ii A, tr B], 

i,j=l i,j=l 



we can compute, based on ( 4.14| ), that 

{trL^-(A),trL'(^)} = tr{L^XA) ? ^'(/^)} 
= tr[r'=''(A,^),Li(A)+L2(^)] =0, k,l>l. 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



This will be used to generate an involutive system of functions defined over the symplectic 
manifold (M^^^jU;^) for any natural number m. 

4.2 An involutive and functionally independent system 

Let us begin to construct an involutive system of polynomial functions by expanding 



det(z^/„ - L(A)) = z.™ - ^z.'"-! + .Ff ^ . . . ^ {-l)"'J^^^\ u = const., (4.18) 
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where 1 < k < m, must read as 



E 

l<Ji<i2<---<ife<m 



L 



^32]! Lj2j2 ' ' ' Lj2jk 



^jkjl ^jkj2 ' ' ' ^jkjk 



, l<k<m. (4.19) 



Here we mention once more that L = {Lij)mxm is assumed. We define bihnear functions Q)^ on 

N 



1 l>0 



Y,{A'^^,B^j)X-'~\ l<i,j< 



where A and B are given by ( |2.40| ), and and are defined as before 

= {4>il,(pi2, - ■ ■ , *i = (V'il, V'i2, • ■ ■ , V'iAf)^) 1 < i < m. 

Then we have 



(4.20) 



(4.21) 



m. 



l>0 

La = Ci + ^{A^^i, B^i)X-^-^ = Ci+ Q^, l<i< 



m. 



l>0 



Therefore, the system of functions J-?'^ is transformed into 



(fc) 



E 



l<Ji<i2<--<ifc<m 



jljl jlj2 jljk 

Cji+ Q X Q X ■■■ Q X 
hh hh j2jk 

Q \ + Q A • • • Q X 



jkji 
Qx 



]kj2 

Qx 



jkjk 

^3k+ Q X 



, l<k<m. (4.22) 



(k) 

A set of more concrete formulas for computing JF^ will be given in Appendix B. Now we further 
expand J^^^ as a power series of 1/A: 

-(fc) _ ^(fc) 



^A = -^1(^1, ■■■,Cm) = Y. Fkiici, • • • , c^)A-', 1 < A; < 



m. 



(4.23) 



l>0 



Based on the formulas of ^1 in Appendix M, it is not difficult to find that 



-^feo — -^fco(ci, ■ ■ ■ ,c„ 



Fki — Fki{ci, ■ ■ ■ , c„ 



E n 



i<ii<i2<---<ife<nip=i 

min(fc,/) fc 

E E E n 

l<il<i2<-<jfc<m r=l l<ii<i2<-<ir<k p=i 

P7^ii,i2,---,ir 
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Pl+P2-\ \-Pr = l~r 

Pl,P2,---,Pr>0 



I > 1, (4.24) 



which are ah polynomials in the canonical variables (pis and tjjis, 1 < i < m, 1 < s < N . 



Theorem 4.2 For all constants ci, C2, • • • , Cm; the polynomial functions in (pis and ipis, 1 < i < 
m, 1 < s < N: Fii{ci, ■ • • ,Cm), 1 < i < m, I > 1, defined by ( ^.24 ), are in involution in pair 



with respect to the Poisson bracket )■ 



Proof: On the one hand, by using Newton's identities on elementary symmetric polynomials 



where 1 < k < m and 



r(2). 



Ci(A) = trL^(A), l<i<m, 



we can have 



-^f = 4'^(Ci(A), C2(A), • • • , a(A)), 1 < A: < m. 
Therefore, we can compute that 

{.^f ,.^^)} = {J-f (Ci(A),C2(A),---,a(A)),.^^)(Ci(/u),C2(^)r--,0(M))} 



(4.25) 



1=1 j 



-^90 (A) 50 (/i) 



The last equality is a consequence of the involutivity of Cj(A), 1 < i < m, shown in ( 4.17 ). On 
the other hand, we have 

{4'\TP}=Y,{Fk^,F,,}X-^^,-K 

i,j>o 

It follows that the polynomial functions Fn = Fii{ci, ■ ■ ■ , Cm) , 1 < i 1^ "rn, I > 1, are in involution 
in pair with respect to the Poisson bracket ([4.5|). I 

Let us now go on to show the functional independence of the polynomial functions Fis{ci, ■ ■ ■ ,Cn 
l<i<m, 1< s < N. 



Theorem 4.3 // all constants ci, C2, • ■ ■ , Cm o,re distinct, then the polynomial functions in (pis 
and ipis: 1 ^ * ^ w,; ^ ^ s < N : Fis{ci, • • ■ , Cm ), 1 < i < m, I < s < N, defined by ^.2^ , are 
functionally independent over a dense open subset o/M^™"^. 



23 



Proof: Let Pq be a point of M^™"^ satisfying 

(pis = e, I < i < m, I < s < N, 



where e is a small constant. Keep ( [4.24| ) in mind, and then at this point Pq, we obviously have 



•rjS2 YJS2 l<ji<j2< - <i,<mg=l P=i 



(4.26) 



l<Jl<i2<---<-'i-l<'" 



where 1 < i, j < m, 1 < si, S2 < N . In the above computation, only the term with r = 1 in the 
expression ( 4.24| ) of Fig contributes to the first-order term of e. Let the matrix B^r be defined 
by 

Qn = (e^f VxTv, Q\P = ArVi, l<hj<N, 
whose determinant is easily found to be 

N 

det(e,v) =n/^i n (^i-^*)- 

i=l l<i<j<N 

Then at the point Pq, the Jacobian of the functions Fj^i with respect to ipjs2 can be computed 
as follows 

d{Fii, • • • , Fin, F21, • • • , F2N, ■ ■ ■ , Fmi, • • • , i^mTv) 



• • ■ ,tplN,tp21, ■ ■ ■ ,'>p2N, ■ ■ ■,fpml, ■ ■ ■ AmN) 



&N ^QGat ^ CiCjQ 

i=2 2<i<j<m 
m 



N 



1=2 
m 



N 



i = l 



1 <z<J <m 



e 



N 



m— 1 



17^2 



m— 1 



i=l 



TV 



CiCjQN 

l<j<jr<m— 1 

det(17„ ® Otv) + 0(e™^+^) 
e'"^(det(0„))^(det(G^))"^ + 0(e™^+^) 

AT 



n 



i<i<j<m i=l l<i<j<N 

where we have used the determinant property of the tensor product of matrices and the deter- 
minant result of the matrix in Appendix This allows to conclude that if the constants 
ci,C2, - • • ,Cm are distinct, the above Jacobian is not zero at Pq when e 7^ is small enough. 
Since the Jacobian is a polynomial function of (pis and ipis, 1 < i < m, 1 < s < N , it is not zero 
over a dense open subset of M?"^^ . Therefore, the functions Fig, 1 < i < m, I < s < N, are 
functionally independent over that dense open subset of M?"^^ . The proof is complete. I 
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4.3 An alternative involutive system to the Fj^'s 





■,Vm) = 1, 








' 'Vjf,, 1 < k < m 




l<jl<i2<-<jfc<m 






■ 5 Vjn) = 0, when k > m + 1 


or k < —1, 



We would like to express the involutive system of the polynomial functions Fig in another way, 
and so we introduce 



(4.27) 



where vi,V2, ■ " are m numbers. Obviously, for m > 2, we have the following relation 

Skivi, - ■ ■ ,Vm) = VmSk-l{vi, - ■ ■ ,Vm-l) + Skivi, - ■ ■ ,Vm-l), k £ Z. (4.28) 

Let us now define 

i-l 

Ell = Fii, Ea = i-iy+'Fii + ^(-iy+h,ici, • • • , c^)^,-^, z, i > 2, / > 1. (4.29) 



From (1^), we can have 



Fa = • • • ,c™)F^-i,z, i,l>l. 

3=0 

Therefore, by Proposition D.2 in Appendix we obtain 



(4.30) 



Eii = Eii{ci, - ■ ■ , c„ 



min(i,Z) 



E E 

l<ii<i2<---<jr<"!. h+h^ ^lj.=l-r 

h,l2,---,lr>0 



M ... J"" 

Si i2 



E 



P1+P2H \-pr=l—r 

pi,P2,---,Pr>0 



{AP^<^j„B^j,) {AP-^<^j„B^j,) ■■■ {APr<^,^,B^,,) 
{AP^^j„B^j,.) {AP^^j„B^j,.) ... (AP'd>,„,5vi/^„) 



(4.31) 



where 1 < i < m and / > 1. Obviously, each En is a linear combination of the -Fjz's, and hence 
{Eiii.,Eji} = holds for all 1 < i, j < m and k,l >1. This means that the polynomial functions 
Eis, 1 < i < m, 1 < s < N , are also in involution in pair. 

In order to show the functional independence of Eis, 1 < i < m, 1 < s < N , similar to the 
proof of Theorem let Pq be a point of M^™^ satisfying (pis = e, 1 < i < m, 1 < s < N, 
where e is a small constant. Then at this point Pq, we have 

dE, 



dip 



-= ed-^Xll-^f,s2 + O(e'), 1 < i,j < m, 1 < si,S2 < N. 



JS2 



(4.32) 
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Hence a direct argument can give rise to 

^(£'11; • • • 1 EiN, E21, ■ • • , E2N1 ■ ■ ■ , Eml, • • • , EmN) 
d{lpu,- ■ • ,'i/'l7V, ■021, • • • ,V'2Af, • • • , V'ml, • ' ' , V'mAr) 

N (4.33) 



i=l l<i<j<N ^<i<j<m 



Therefore, if ci, C2, • • • , Cm are distinct, the above Jacobian is not zero at Pq when e 7^ is smah 
enough. This imphes that the functions Eis, 1 < i < m, 1 < s < N , are functionally independent 
over a dense open subset of M^*"^. 

Let us sum up these results in the following theorem. 



Theorem 4.4 AH polynomial functions in (pis and ipis, 1 < i < m, l<s< N : Eii{ci, • • • , Cmj, 
1 < i < m, I > 1, defined by ( 4-3i ), are in involution in pair with respect to the Poisson 



bracket ^.^j ) for all constants ci,C2, ■ ■ ■ ,Cm- Moreover, among them the polynomial functions 
Eis{ci, ■ ■ ■ ,Cm), 1 < i < m, 1 < s < N , are functionally independent over a dense open subset 
of R^"^^ for distinct constants ci , C2 , • • • , Cm • 

Note that all polynomial functions Fn are also linear combinations of the En^s. The above 
theorem actually shows us an alternative to the involutive and functionally independent system 
of the polynomial functions Fig, 1 < i < m, 1 < s < N . The Sj^'s have the compact form for 
the constants ci, C2, • • • , Cm, and thus it is more convenient to deal with them. 

5 Liouville integrability and involutive solutions 

Let us now turn to establish the Liouville integrability of the obtained constrained flows, and 
to present involutive solutions of the TV-wave interaction equations in both 1 + 1 and 2 + 1 
dimensions. The involutive system of the polynomial functions 

Fis = Fis{ci, ■■■ , Cm), 1 <i <m, 1 < s < N, 

alternatively 

Eis = Eis{ci, - ■ ■ ,Cm), l<i<m, l<s<N, 
will play an extremely important role in the following discussion. 

5.1 Liouville integrability of the constrained flows 

For the 1+1 dimensional case, we have the matrix Lax operator as defined by ( p. 49 ) and ( 2.501 ), 

(A) = (A; 71 , • • • , 7n) = Ci (71 , • • • , 7n) + i?i (A) , 



I.e. 
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where Ci and Di{X) are given by (|2^ ). Note that 



H ^Ij, 'i- j < n. 



According to Theorem 4.2 and Theorem |4.3| for the case m = n and Q = 7^, 1 < i < n, we know 
that -Fis(7i, • ■ ■ ) 7n), l<^<'ra, 1 < s < N, defined by ( 4.24| ), are functionally independent 
over a dense open subset of M^"^ and in involution in pair with respect to the Poisson bracket 
O), i.e., 



i=l 



oofm,2nN\ 



Theorem 5.1 Let 71,72, • • • ,7n n distinct numbers. Then the spatial constrained flow (2.4j) 
and the temporal constrained flow ( 2.44 ) of the 1 + 1 dimensional J\f -wave interaction equations 
are Liouville integrable Hamiltonian systems, which possess involutive and functionally 
independent integrals of motion 

Fis{-fi,---,-fn), l<i<n, l<s<N, 



defined by ( 4-24 ) in the 



case 



m = n, Ci = ji, 1 < i < n. 



Proof: From the necessary Lax representations of the spatial constrained flow ( |2.43| ) and the 
temporal constrained flow ( ^.44 ): 



(L«(A)), = [[/(n,A),LW(A)], (L«(A))i, = A), L«(A)], 

which are shown in Theorem 2.1, we can obtain |26| 

(L«(A)r),, = [U{u,X),{L^'Hm, {L^'Hmt, = [V^'\u,X),{L^'HX)n iJ > 1, 

and thus we have 

{tv{LW{X))% = ti{{LW{X))% = tr[[/(n. A), (L«(A))i = 0, i > 1, 
(tr(L(i)(A)y)t, =tr((L(i)(A)y)i, =tr[y(i)(ti,A),(L(i)(A)y] =0, j > 1. 



Therefore, .^{^^^(71, • • • , 7n) are all generating functions of integrals of motion of ( 2.43| ) and 
( |2.44| ) in the light of the expression ( |4.25[ ) determined by Newton's identities. It follows that 
F'isilir ' ' tIu), ^^il^n, 1<s<N, are all integrals of motion of the spatial constrained 
flow ( 2.43| ) and the temporal constrained flow ( 2.44| ) . Note that all constants 71 , 72 , • • • , 7n are 
distinct. Therefore, Theorem |4.2| and Theorem |4.3| in the case of m = n and q = 7^, 1 < i < n, 
together with Theorem 2.1, show that the spatial constrained flow ( 2.43| ) and the temporal 
constrained flow ( p. 44 ) are Liouville integrable Hamiltonian systems, which posses the involutive 
and functionally independent integrals of motion Fisi^ji, ■ ■ ■ ,jn), ^ < i 1^ n, 1 < s < N. The 
proof is finished. I 
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We remark that from the Lax representations shown in Theorem 



2.1. 



, we have 



K„-L«(A)),. = [[/(n,A),i//„-L«(A)], 
[vl^ - L«(A))t, = [y«(n, A), uln - L«(A)] 

for any constant v. It fohows H] that det(i//„ - L(i)(A) is a common generating function of 
integrals of motion of the constrained flows (|2.43|) and ( ^.44|) , and thus so are ^^'^^(71, • • • ,7n), 
1 < k < n. This is an alternative proof for showing that ^^'^^(71, • • • ,7n), 1 < k < n, are the 
generating functions of integrals of motion of ( 2.43| ) and ( ^.441) . 



For the 2 + 1 dimensional case, a completely similar argument can give rise to the following 
theorem on the Liouville integrability of the constrained flows ( p.23[ ), ( |3.24 ) and ( p.25 ) of the 
2 + 1 dimensional A/"- wave interaction equations (p.4[). 



Theorem 5.2 Let 5i, • • • be n + 1 distinct numbers. Then all three constrained flows 

( 3.2^) , ( 3.2^ ) and ( 3. 2^ ) of the 2 + 1 dimensional M -wave interaction equations (3^) are Liou- 
ville integrable Hamiltonian systems, which possess the involutive and functionally independent 
integrals of motion 

Fis{Si, ■■■ ,6n, Sn+i), l<i<n + l, l<s<iV, 
defined by fi4-24) in the case 



m 



n + 1, Ci = 6i, 1 < i < n + 1. 



5.2 Involutive solutions of the A/^-wave interaction equations 

We would like to show that the constrained flows provide involutive solutions to the M-wave 
interaction equations in both 1 + 1 and 2 + 1 dimensions. For the 1 + 1 dimensional case, we 
have the following result. 



Theorem 5.3 If (j)is{x,ti) and ipis{x,ti), l<i<n, l<s<N, solve the spatial constrained 
z tempora 

Uij{x,ti) 



flow ( ^-4^ 0-nd the temporal constrained flow ( ^-441 ) simultaneously, then 

ai — a 



-{^i{x,ti),B^j{x,ti)), l<i^j<n, 



(5.1) 



7i - 7i 

with ^i{x,ti) and ^i{x,ti) being given by 

$i(x,ti) = ((/>ii(x,ti), • • • ,(/)iAr(x,ti))^, '^i{x,ti) = • • • ,-i/'i7v(x,ii))^, 1 < i < n, 

solve the 1 + 1 dimensional J\f -wave interaction equations ( 2.2^ . 



Proof: Note that the 1 + 1 dimensional TV- wave interaction equations ( 2.22| ) is the com- 
patability condition of the spectral problem ( |2.1| ) and the associated spectral problem ( p. 17 ) 



with m = 1 or the adjoint spectral problem ( 2.27 ) and the adjoint associated spectral problem 
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( 2.2§| ) with m = 1 for whatever potential u. Therefore, the 1 + 1 dimensional TV- wave interaction 
equations ( ^.22 ) is also the compatability condition of the spatial constrained flow ( |2.43| ) and 
the temporal constrained flow ( 2.44D under the constraint ( |2.41 ). Now (/)js(rc,ti) and 
l<i<n,l<s<A^, are assumed to solve ( |2.43 ) and ( |2.44| ) simultaneously, and thus the 
potential defined by ( ^.iD must satisfy the compatability condition of the spatial constrained 
flow (|2.43| ) and the temporal constrained flow ( ^.44 ). This means that the potential defined by 
(|5.1| ) must be a solution to the 1 + 1 dimensional A/"- wave interaction equations ( 2.22| ). The proof 
is finished. I 

We remark that a direct computation can also show the above theorem. For the 2 + 1 
dimensional case, a similar deduction can give rise to the following theorem. 



Theorem 5.4 If 4>is{x,t) and ipis{x,t), 1 < i < n + 1, 1 < s < N, solve the constrained flows 
^3.2$) , ( ^.24 ) and (\3.23^ simultaneously, then 



Pij{x,y,t) 



Ji J' 



5i - 6j 



^{^i{x,y,t),B^j{x,y,t)), I < i ^ j < n, 



Ki — K 

qij{x,y,t) = -j—j^{^i{x,y,t),B-^j{x,y,t)), 1 < i ^ j < n, 



(5.2) 



with $j(x,t) and '^i{x,t) being given by 

$i(x,t) = {(j)ii{x,t), ■ ■ ■ ,(t)iN{x,t)f , '^iix,t) = {ipii{x,t), ■ ■ ■ ,tl)iN{x,t)f , I <i <n+l, 
solve the 2 + 1 dimensional M -wave interaction equations ( ^■4 )- 



Also, one can find that 



^ {<l>i,B^n+l), 9i = T \ {^n+l,B^i), l<i<n 



(5.3) 

(>i — On+1 Oi — On+1 

provide a solution to the Lax system (|3.lD and the adjoint Lax system ( |3.5D with the potentials 
given by (|5.2| ). What's more, ( ^.2| ) and (|5.3| ) automatically satisfy our first symmetry constraint 



In the following theorem, the solutions given in Theorem 5^ and Theorem 5^ are shown to 
be involutive. 



Theorem 5.5 The Hamiltonians Hf and H^^ of the constrained flows in 1 + 1 dimensions, 
defined by ( 2. 55 ) and ( \2.54 ), are the second-order polynomial functions of Eii{'yi, ■ ■ ■ ,jn), 1 < 
i <n, I = 1,2, and thus they commute, i.e., 

{Hf,H\^] = d, (5.4) 

where the Poisson bracket {•, •} is defined by ^2.4V - The Hamiltonians H2, and H2 of the 
constrained flows in 2 + 1 dimensions, defined by ( ^. 3(\ ), ( 3.31 ) and ( 3. 3^ ), are also the second- 
order polynomial functions ofEii{6i, ■ ■ ■ ,5n,Sn+i), 1 <i < n+1, I = 1,2, and thus they commute 
with each other, i.e., 

{m,Hl} = {HI, Hi} = {Hi Hi} = 0, (5.5) 



29 



where the Poisson bracket {•, •} is defined by ( 3. 28 ). 



Proof: Directly from the explicit expression (|43ll ) of the Eis's, we have 



E^l = J2c'-\<^j,B^j), l<i<m, (5.6) 



E^2 = X]c}-l(^<I>„i?^',) 



l<j<k<m 



m J-1 



E^r'^i- 7^(^j,^*j)(^fc,^*fc>, l<i<m, (5.7) 

hk=l 



j=l j,k=l '^^ 



where the fj's are defined as follows 



rn ^ 

£j = {A^j,B^j) + ^ -^—{q>.,B^k){<!>k, B^j), l<j<m. (5.8) 
fc=i ^'^ 

Now solving ( |5.6D for (<I>j, 1 < i < m, leads to 

= (n ) J]](-l)"'"-''sm-i(ci,--- ,Ci-i,Ci,Ci+i,--- ,c,„)Sji, 1 < i < m, (5.9) 



where the Sj's are defined by ( 4.27] ) and q means that q does not appear. Therefore, each 



{^i,B"^i) can be expressed as a linear combination of En, 1 < i < m. Similarly, solving (|5 
for Sj, 1 < j < m, leads to 

m ^ m 

Si = (TT ) y'(-l)™"-^Sm-j(ci,---,Ci_l,Q,Ci+l,---,Cm) X 

r=l j=l 

m J-1 



Ej2+ V — ^ l<i<m. (5.10) 

:,l=l 
k^l 



k,l=l 



This expression together with ( |5.9| ) implies that each £j can be expressed as a linear combination 
of £"4! and Ei2, 1 < i < m. 

In the 1 + 1 dimensional case, we have m = n, cj = 7j, 1 < i < ?^. Hence 

n ^ 

£j = {A^j^B^j) + -^—{^j,B^k){^k,B^j), l<j<n. (5.11) 

k=i 
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The Hamiltonians Hf and H^^ in Theorem 2T can be easily expressed as 

n n 



(5.12) 



k=l k=l 

where the £kS are defined by ( |5.11 ). 

Likewise, in the 2 + 1 dimensional case, we have m = n + 1, Cj = 5j, 1 < j < n + 1. Hence 



n 



6i — 6k 
k=i ^ ^ 



+ 



6j — Sn+1 



{<^j,B^n+l){'^n+l,B^j), l<j<n, 



£n+l = {A^n+l,B^n+l) + Yl ' 



k=l 



>n+l 



{^n+l,B-^k){^k^B-^n+l) 



(5.13) 
(5.14) 



The Hamiltonians -fff , i^l and H2 in Theorem 3.1 can be expressed as 

n n n 



(5.15) 



k=l k=l k=l 

where the S^s are defined by ( |5.13 ). 

Therefore, Hf and are linear combinations of £'^^(71, • • • , 7„), 1 < i < n, Z = 1, 2, and 
H2, H2 and H2 are linear combinations of En^Si, ■ ■ ■ , 5n, 6n+i), 1 < i < n + 1, Z = 1, 2. It follows 
from Theorem iA that Hf and are in involution, and i/f , H2 and H2 are in involution in 
pair, too. The proof is finished. I 

We remark that a direct computation can also give a proof for the involutive property of 
the Hamiltonians of the constrained flows in both 1 + 1 and 2 + 1 dimensions. Only a new set 
of equalities 

h cycle(z, j,k) =0, l<i,j,k <n, 

Cj Ci Ck Ci Ck Ci Cj Ci 

has to be utilized, where Cj, 6j, and Ci, 1 < i < n, are arbitrary constants. This just needs a 



direct check, too. However, the proof of Theorem 5.5 also gives rise to the explicit expressions 
for all Hamiltonians of the constrained flows in both 1 + 1 and 2 + 1 dimensions, in terms of the 
integrals of motion Eis . 

Now if we denote the Hamiltonian flows of the spatial constrained flow ( p. 43 ) and the 



temporal constrained flow ( 2.44 ) by Qx ^ and ^ respectively, then the above theorems present 
a kind of involutive solutions to the 1 + 1 dimensional A/"- wave interaction equations ( p. 22 ): 



Uij{x,tl) 



ai — aj . HI h\ 



jO 

li - Ij 

^-{Qt ' gx ' $io, gt " gx ' B^jo), l<t^j<n, 

li - Ij 



(5.16) 
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where the initial values and of <I>j and can be taken to be any arbitrary constant vectors 
of the Euclidean space . Similarly, if we denote the Hamiltonian flows of the constrained flows 
(|3.23| ), (|3.24|) and (|3.25D Qx^ , Qy^ and ^ respectively, then the above theorems present a 



kind of involutive solutions to the 2 + 1 dimensional A/"- wave interaction equations (|3.4| ) : 

Pij[x,t) = — -^{g^^gy^g^^'^io^g^r^gy^gt B^jo) 

- - 

y^9t ^9x^'^io,gy^gt ^gx^B^jo) 





- ^ 


Ji- 




Si- 


-s - 


Ji- 


-Jj, 


Si- 






1 < I 


Ki 




Si 








Si 









t jjx jjy _ jq^x jjv _ 

'gx^gy^^io,gt ^gx^gy^B^jo) 



-L^'i i ' H^ H^ H^ — H^ H^ — 

qij{x,t) = ^ — Y^{gx^ gy'' gt^^io,gx'' gy^ gt'^B^jo) 

ZjV Ljt TJX _ UV Zjt TJX 

J^O T. J^O -f^o 



(5.17) 



'■gt ^gx^^io,gy^gt ^gx^B^jo) 



H^ H^ H^ ~ H^ H^ H^ ~ 

h ^gx ^ gy ' ^io, gt ^9x ' gy b^jq) 



Si - Sj 

= l<i^j<n, (5.18) 

where the initial values and ^io of and ^i can also be taken to be any arbitrary constant 
vectors of the Euclidean space M^. 

Note that all constrained flows in both 1 + 1 and 2 + 1 dimensions are Liouville integrable, 
and that the initial values of <l>j and ^i, 1 < i < n, can be arbitrarily chosen. Therefore, 



together with Theorem 5.1 and Theorem the above involutive solutions also show us the 
richness of solutions and the integrability by quadratures for the A^-wave interaction equations 
in both 1 + 1 and 2 + 1 dimensions. Of importance is of course that binary symmetry constraints 
decompose the A/'-wave interaction equations in both 1 + 1 and 2 + 1 dimensions into finite- 
dimensional Liouville integrable Hamiltonian systems, and the resulting involutive solutions 
present the Backlund transformations between the TV-wave interaction equations in both 1 + 1 
and 2 + 1 dimensions and these finite-dimensional Liouville integrable Hamiltonian systems. 



6 Conclusions and remarks 

We have introduced a class of special symmetry constraints, ( |2.38D in the 1 + 1 dimensional case, 
and ( 3.18| ) and ( p. 201 ) in the 2 + 1 dimensional case, for the A/'-wave interaction equations in 



both 1 + 1 and 2 + 1 dimensions. These symmetry constraints nonlinearize the n x n spectral 
problem and adjoint spectral problem, ( 2.34| ) and ( 2.35| ), and the (n + 1) x (n + 1) spectral 



problem and adjoint spectral problem, ( ^.11 ) and ( 3.12| ), into finite-dimensional Liouville inte- 



grable Hamiltonian systems, and decompose the AA-wave interaction equations in both 1 + 1 and 
2 + 1 dimensions into these finite-dimensional Liouville integrable Hamiltonian systems. A gen- 
eral involutive and functionally independent system of the polynomial functions Fis(ci, • • • , Cm), 
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1 < i < m, 1 < s < A^, or alternatively Eis{ci, ■ ■ ■ , Cm), 1 < i < m, 1 < s < N , associated with 
an arbitrarily higher-order matrix Lax operator, was presented and used to show the Liouville 
integrability of the resulting constrained flows. The nonlinear constraints on the potentials, re- 
sulting form the symmetry constraints, also provide us with a class of Backlund transformations 
from the AA-wave interaction equations to the obtained flnite-dimensional Liouville integrable 
systems. The involutive solutions to the A/'-wave interaction equations are given through the 
constrained flows, and thus the integrability by quadratures has been exhibited for the A^-wave 
interaction equations. The special case with F = Wq, i.e, diag(7i, • • • , 7„) = diag(/3i, • • • , of 
two reductions of n = 3 and n = 4 in 1 + 1 dimensions presents all results established in |3l|, ^ . 

We point out that for a more general matrix Lax operator L = C + D with any constant 
matrix C = {cij)mxm and the matrix D defined by ( |4.2| ), the r-matrix formulation ( 4.13| ) still 
holds. Therefore, an involutive system of polynomial functions can be generated, but we do not 
know what conditions on the matrix C can ensure the functional independence of that involutive 
system. We are also curious about other examples of higher-order matrix Lax operators which 
lead to involutive and functionally independent systems. Our crucial techniques to present the 
involutive and functionally independent system Fis, 1 < i < m, 1 < s < N, are the r-matrix 
formulation, Newton's identities on elementary symmetric polynomials, and the determinant 
property of tensor products of matrices; and the whole process of their applications provides an 
efficient way to show the involutive property and the functional independence. 

Of course, one of the important results in binary nonlinearization is the integrability of soli- 
ton equations by quadratures, which implies that one can integrate soliton equations themselves 
by quadratures. However, the potentials obtained by symmetry constraints can be proved to 
belong to a kind of finite-gap type solutions containing multi-soliton solutions, and thus they 
may not present solutions to given initial value and/or boundary problems of soliton equations. 
It is a challenging problem to establish a general theory of complete integrability for nonlinear 
differential and differential-difference equations, which should state what mathematical proper- 
ties the equations must possess so that their solutions to initial value and/or boundary problems 
can also be determined by quadratures. 

Symmetry constraints yield nonlinear constraints on potentials of soliton equations, and put 
linear spectral problems (linear with respect to eigenfunctions) into nonlinear constrained flows 
(nonlinear again with respect to eigenfunctions), which makes it more complicated to solve 
soliton equations. However, since spectral problems are overdetermined, one needs additional 
conditions (compatability conditions) to guarantee the existence of eigenfunctions of spectral 
problems. The symmetry property brings us the Liouville integrability for nonlinear constrained 
flows. Thus, symmetry constraints make up for the disadvantage of nonlinearization in manip- 
ulating binary nonlinearization. Of special interest in the study of symmetry constraints are 
to create new classical integrable systems [^], which supplement the known class of integrable 
systems |56], and to expose the integrability by quadratures for soliton equations by using con- 
strained flows [33|. 
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The idea of binary nonlinearization is quite similar to that of using adjoint symmetries to 
generate conservation laws for differential equations, both Lagrangian and non-Lagrangian 
In binary nonlinearization, we adopt adjoint spectral problems to formulate Hamiltonian struc- 
tures for constrained flows so that finite-dimensional Liouville integrable systems result. Note 
that there exist also some special symmetry constraints which do not yield Hamiltonian struc- 
tures with constant coefficient symplectic forms, including both canonical and nono-canonical 
ones, for constrained flows |4^. Therefore, it will be particularly interesting and important to 
classify symmetry constraints which exhibit Hamiltonian structures with constant and variable 
coefficient symplectic forms for constrained flows. 
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A Non-Lie symmetries 



Proposition A.l If (1)^"' and tp^"' , 1 < s < N, satisfy ( ^.34\ ) and i2.33i ), then the vector field 



is a symmetry of the 1 + 1 dimensional -wave tnteractton equations 



Proof: It is required to show that 

JV N 

(SP, 5Q) = {[Uo, ^ M.^f^V^^^^l, [Wo, J2 M.-^'^V'^^^]) 



satisfies the hnearized system (2.25). By using (2.34) and (2.35), we can first compute that 

N JV JV 



3 = 1 



3=1 3=1 

JV AT 

s = l s = l 

JV 

= E^4V^<^'KA.),0(='^(=n 

s = l 

JV JV 



and similarly, we can have 



s = l s = l 

Thus, noting the Jacobi identity, it follows that 

JV 



{SP)t,-{SQ), = [[/o,(Em30'^V'^'^)^i]-[M^o,(Em.0'^'V''^'^W 

s = l s = l 

JV 

= E A«M.([f/o, [Wo, -^t^V'^'"^] - [W^o, [;7o, 0^^'^'^]) 

3 = 1 

JV JV 

+[f/o,[m,E/^-'^''V'^'^] - [m^o,[;7i,Em.0^^V^^^''] 

3 = 1 3 = 1 

JV JV 

= [Uo, [m, E M3</>'^V'^'^] - [Wo, [(7i,E M30^^V^^'^], 



(A.l) 



(A.2) 



where 6P and are defined by (A.2). Then, again noting the Jacobi identity, we can have 

{SP)t, - (5Q). + [Ui , 5Q] + [5P, W^] 

JV JV 

= E M3<^'''^'''^, [Uo, Wi]] - E M30<' V'^, [Wo, Ui]] = 0, 

3 = 1 3 = 1 

in the last step of which we have used [Uo, Wi] = [Wo, Ui]. The proof is finished. I 

All of the symmetries presented in this proposition are not Lie point, contact or Backlund symmetries, since 
they can not be written in terms of the potentials Uij and their spatial derivatives. 
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B Formulas for computing JT! 



(k) 



Immediately from the expressions of J-^^'' in (4.22), we can obtain the following more concrete formulas for 
computing T^^: 
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C The determinant of il. 



The following proposition has been used while showing the functional independence of the polynomial functions 
Jis(ci, • • ■ , Cm), 1 < i < m, 1 < s < N, which is of interest itself. 
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Proposition C.l Let m > 2, and ci,C2, ■ ■ ■ ,Cm be constants. Then 



det(Sl„ 



m 


2<i<j<m 


2<z<j<fc<m 


m 


rri 

i=l 


l<i<j<m 
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m 

i=i 



1 E 



E 



CiCjCk 



i — 1 l<i<j<m — l l<i<j<k<m~l 



n 



n 



(C.l) 



l<i<j <m 



Proof: We prove this pr opos ition by tlie principle of mathematical in duc tion. It is obvious that (C.l) is true 
when m — 2. Suppose that (C.l) is true when m — I. Let us verify that ( |g.l| ) is also true when m — I + 1. Note 
that 



E 

(Ci ~Cj) 



l<n<»2<--<ifc<' + l l<^^<^2<■■■<^k<' + ^ 



E 



E 



CiiCi2 ■ • ■ Ci^_i, l<i,j<^ + l, 1 < fc < 



I<il<i2<-<ifc_l<! + l 
il.»2,---,ifc-l^i,j 



For each 2 < j < Z + 1, we subtract 



CijCi2 ■ ■ • Ci^^i X the first column of det(f2i+i) 

2<ii<i2<---<ij-l<'+l 



from the jth column of det(r2i+i), and then we have 
det(r2i+i) 
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in the last step of which we have used the inductive assumpt ion. This means that ( C.l ) is also true when m — l + l, 
i.e., the inductive step is satisfied. Therefore, the formula (CT) is always true by the principle of mathematical 
induction. The proof is finished. I 



D Two identities on symmetric polynomials 



Let the Sj's be symmetric polynomials defined by (4.27) 



Proposition D.l For any integers r and i with i > r > 1, and any numbers ci, ■ ■ ■ , Cr, we have 

1, if i = r, 
0, if i > r. 



Y^{-lYs,{Cl,---,Cr) 



j=0 



/iH \-lr—i — r — j 

ll,---,lr>0 



(D.l) 



Proof: Use the principle of mathematical induction on r. Whe n r — 1 and i — 1, the left-hand side of (D.l) 
is 1. When r — 1 and i > 1, the left-hand side of (D.l) is 0. Hence (D.l) holds when r — 1. 

Now suppose that ( |D.i| ) holds when r = k, i.e., 

^-^ I' 1, ifj = fc, 



ilH \-l^.—i-k-j 



^(-iy.,(ci,...,c,) 

Then, when r = k + 1, the left-hand side of ( |D.1[ ) is 

i-k-l 



0, if j > k. 



j=0 



'iH hife + l=i-fc-j-l 

il,---,'fc + l>0 



By using (4.26), it equals to 

i—k—l i—k—j—1 



E (-1)' E ci^:i'''sj-^{c„---,ck) E 



^'l Jk 
1-1 L-. 



j=0 ir+l=0 
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il....,!fc>0 



+ E (-1)' E c[':isAci,---,ck) 

3 = ir + l=0 



E 
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E 4Ti E (-1)^^^^^.(ci,---,c. 

i,.+ l=0 j=0 



E 



^1 ^fc 

'iH h!fc=»-*:-j-'fc + i-2 



+ E 47/ E E 

ir+l=0 3=0 



iiH hifc=*-fe-i-ifc+i-i 

il,---,!fc>0 



where an empty sum is understood to be zero. 



(D.2) 



(D.3) 



(D.4) 



When i = fc -I- 1, it is easy to see that (D.4) equals to 1. If i > fc -f 1, then by (D.2), the first sum equals to 



and again, by (D.2), the second sum equals to 

Jk+i I 
fc+i 



+ 1 =i-k-l 



'-k + 1 I 



-k + l 



(D.5) 



(D.6) 



Hence (D.4) equals to if i > fc -f 1, which implies that (D.l) holds when r — k + l. Therefore, (D.l) always holds 
by the principle of mathematical induction. The proposition is proved. I 
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Proposition D.2 For any integers m, r, i with i>r + l>2, m numbers ci, ■ ■ ■ , Cm, and r integers ji, • ■ ■ ,jr 
with 1 < ji < ■ ■ ■ < jr < rn, we have 



^(-1)' ^Sj(ci,---,c„) 

j=0 



c ■ ■ ■ c = 



ylr — i-r-j 



l<pl<--< p^_^<m 

Pa7^j/3 for all a, (5 



(D.7) 



Proof: Without loss of general ity, s uppose that ji — i when i = 1, ■ ■ ■ ,r, since each Sj (ci, • • ■ , c,„) is symmetric 
with respect to ci, ■ ■ ■ ,Cm- Then, (|D.7|) becomes 



^(-1)' •'sj(ci,---,c™) ^ 
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(D.I 



j=0 



ilH \-lr — i — r — j 

ll,---,lr>0 



r + l<pi<---<p^_^<m 



Obviously, for any fixed j with r + 1 < j < m, both sides of (D.8) are linear with respect to Cj 



We use the principle of mathematical induction on i to prove (D^). When i — r + 1, both sides of ( D.8 ) 
equal to Cr+i H he™. 



Suppose that (D.8) holds when i = k {k > r). Then, when i = k + 1, the left-hand side of (D.i) reads as 



R := 'Y.^i-ir+'-^-^s,ic,,---,< 



E 



J=0 
fc-r 



il+'-' + ir^fc + l-r-j 



il+--- + 'r = fc-»--i 

ll,---,ir>0 



Then by (4.28), we have 



j=0 



il + '" + ir = fc-r-j 



By the inductive assumption, it becomes 

m 

dc 

Hence we obtain 



OR _ 



r + l<pi<---<p^_j.<m — l 



R— ^ Cpj • ■ ■ Cpj._,,Cm + 7?l(ci, • ■ ■ , Cm_l), 

r + l<pi<---<pf._^.<m — l 

where Ri is a polynomial. Since R is symmetric with respect to Cr+i, ■ ■ ■ , Cm, we have 
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+ -Ro(ci, ■ • ■ ,Cr. 



'■ + l<Pl<---<Pfc + l-r<»" 



where by setting Cr+i = ■ ■ ■ = Cm = in (D.9), Ro is determined to be 



k + l-r 



Cl 



1 5 ^ry 



E 



3=0 



By Proposition D.l, i?o =0 since k + 1 — i > r. Hence 

R = 
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ilH f-!r = fe + l-r-j 
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(D.9) 

(D.IO) 

(D.ll) 
(D.12) 

(D.13) 
(D.14) 
(D.15) 



which implies that (D.8) holds when i = fc + 1. Therefore, (D.8) holds for alH > r by the principle of mathematical 
induction. The proof is completed. I 



The identity (D.7) is needed in presenting an alternative involutive system Eis's to the Fia's in the subsec- 



tion 4.3, 
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